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ABSTRACT 



Aims. High-energy irradiation of the circumstellar material might impact the structure and the composition of a protoplanetary disk 
and hence the process of planet formation. In this paper, we present a study on the possible influence of the stellar irradiation, indicated 
by X-ray emission, on the crystalline structure of the circumstellar dust. 

Methods. The dust crystallinity is measured for 42 class II T Tauri stars in the Taurus star-forming region using a decomposition 
fit of the 10 yum silicate feature, measured with the Spitzer IRS instrument. Since the sample includes objects with disks of various 
evolutionary stages, we further confine the target selection, using the age of the objects as a selection parameter. 
Results. We correlate the X-ray luminosity and the X-ray hardness of the central object with the crystalline mass fraction of the 
circumstellar dust and find a significant anti-correlation for 20 objects within an age range of approx. 1 to 4.5 Myr We postulate that 
X-rays represent the stellar activity and consequently the energetic ions of the stellar winds which interact with the circumstellar disk. 
We show that the fluxes around 1 AU and ion energies of the present solar wind are sufficient to amorphize the upper layer of dust 
grains very efficiently, leading to an observable reduction of the crystalline mass fraction of the circumstellar, sub-micron sized dust. 
This effect could also erase other relations between crystallinity and disk/star parameters such as age or spectral type. 



Key words, circumstellar matter - 
stars 



1. Introduction 



stars: pre-main sequence - stars: formation - planetary systems: protoplanetary disks - X-rays: 



The evolution of the dust in a protoplanetary disk is one of the 
key subjects in the overall research on mechanisms of planet 
formation. As we now know, dust in young circumstellar disks 
differs significantly from the dust in the interstellar medium 
(ISM). There is evidence for grain growth from the typical ISM 
and sedimentation in the vertical directio n (or dust settling) of 
a disk (see various references such as Ro dmann et al.l l2006l 
ISiciha-Aguilar et alj IIoOTI iFurlan et a l. 2006). While the dust 
grains in the ISM and in molecular clouds are amorphous, pro- 
toplaneta ry disks contain an increased content of crystalline sil- 
icates ( Bouw man et al. ' 2008). This transition of the dust grain 
structure during the star-forming process is poorly understood 
and could be important in the later scenario of planet formation. 

However, as many authors have pointed out (e.g. 
Watson et al][2 009. Sicilia-Aguilar et al. 2007, van Boekel et all 



2005L ISchegerer et alj 12006 ). no definitive connection have 



been found so far between the properties of the disk or of the 
central object and the crystalline mass fraction of the dust disk. 
From the point of view of standard dust processing scenarios 
for protoplanetary disks, this conclusion is surprising. We 
expect the crystallization of the dust grains to occur due to 



Send offprint requests to: 
glauserOastro . phys . ethz . ch 



M. 



Glauser, e-mail: 



thermal annealing or evaporation and recondensation processes 
either close to the star or within accretion shocks (HenniiiJ 
2008). Radial mixing may transport t he crystalline grains to 
more distant regions (e.g.. !Gail[2004h . Therefore, we expect 
an evolutionary trend for the crystalline mass fraction and/or 
correlations with stellar parameters such as the bolometric 
luminosity, the photospheric temperature, the accretion rate or 
the disk/star mass ratio. The fact that no such relation has been 
found raises the question on alternative me chanisms , contro lling 
the process o f crystallization. Kessler-Sil acci et al] ( l2006h and 
I Watson et alj ('2009) suggested that the crystallizing process 
might be dominated by the impact of X-ray irradiation which 
destroys the crystalline structure of the dust grains. 

Young stars are very strong sources of X-rays. A typical 
T Tauri star emits between 10^'^ and 10"" erg s"' in the soft 
(0.1 - 10 keV) X -ray band, i.e ., 2-4 orders of magnitude more 
than the Sun (seel Gude]|2004l for a review of stellar X-ray radia- 
tion). The radiation is thought to be mostly coronal, originating 
from hot (1-20 million K), magnetically trapped plasma above 
the stellar photosphere, in analogy to the solar coronal X-ray ra- 
diation. 

There is little doubt that X-rays have some impact on the 
gas and dust in circumstellar disks, at least relatively cl ose to 
the sta r and at the disk s urface . Fo r example, |lg ea & Gl assgoldl 
(Il999h . iGlassgold et all (|2004|) or lErcolano et al., (.2008) com- 
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puted detailed models for radial distances between 0.5-10 AU 
that indicate efficient ionization of circumstellar disks by X- 
rays and also heating of the gaseous surface layers to several 
thousands of K elvin. Complicated ch e mical networks are a con - 
sequence (e.g.. ISemenov et aLll2004l or lllgner & Nelsonll2006h . 
A direct evidence for these processes is suggested from the 
presence of strong line radiation of [Ne n] at 12.8 //m detected 
by Spitzer in many T Tauri st ars (e.g., Pascucci et al.l 120071 
iHerczeg et al.ll200"7l or lLahuis et al. 2007) which in some cases 
may be triggered by shocks dvan Boekel et alj2009l) . This transi- 
tion req uires ionization and he ating of the ambient gas to several 
1000 K (iGlassgold et al ."2007^. 

Magnetic energy release events, so-called flares, occurring 
in the same stellar coronae can increase the X-ray output up to 
hundreds of times, but as we know from solar observations, such 
events are also accompanied by high-energy electr ons, protons 
and ions ejected from the Sun. lFeigelson et al] 12002) speculated 
that the expected elevated proton flux around T Tauri stars leads 
to isotopic anomalies in solids in the accretion disk, as suggested 
from measurements of meteoritic composition for our early solar 
system (Caff'e et al. 1987). 

The destructive impact of high-energy irradiation on crys- 
talline structures by io ns has been demo n strated in laborator y 
measurement s by, e .g., Jager et akl ( 2003b . iBringa et"aLl (|2007|) . 
iDemvk et all (1200 Ih ICarrez et al.1 (12002'). mainly for low ener- 
getic cosmic rays (E ^50 keV) but rarely also for lower energies 
typical of stellar winds. 

In this study we look in particular for correlations between 
the crystalline mass fraction and stellar properties related to X- 
ray emission by deriving these parameters for T Tauri stars in the 
Taurus-Auriga star formation region. We present in Sect. |2] the 
target sample and some aspects of the data reduction, describe 
in Sect. [3] the methodology of measuring the crystalline mass 
fraction based on decomposition fits to the lOfim silicate feature, 
present the derived values in Sect. |4] and bring them into context 
with X-ray parameters in Sect.|5] Our conclusions are presented 
in Sect.m 



2. Data sample and data reduction 

We use the sample of objects in common to two recent surveys of 
the Taurus-Auriga star-forming region. The first survey was ob- 
tained by the Spitze r InfraRe d Spectrograph (IRS) and was pub- 
lished previously bylFurlan"et al. (2006) and further analyzed by 
I Watson et al.l (l2009l) . The second survey was perf ormed in the 
X-ray range with XMM-Newton as described by iGiidel et"an 
(l2007h . While the former survey provides information on the 
dust properties of the circumstellar disks, the latter allows the 
investigation of stellar X-rays. We focus only on the disk sur- 
rounded (Class II, as listed in iGiidel et al J 120071) T Tauri stars 
which appear in both surveys and show significant emission in 
the 10 fim silicate feature. Table [T] provides an overview of the 
objects used for this study and their properties derived in the 
framework of the XMM-Newton survey. 

The Spitzer IRS data were obtained during the observing 
campaigns 3, 4 and 12 using mainly the low resolution channel. 
Two exposures per object were obtained in different nod posi- 
tions allowing the subtraction of the background by subtracting 
the two spectra, Fy^i{A) and Fy^2{^), from each other. The fluxes 
are then averaged over the two observations and the difference 
of the two spectra d{A) = Fy i{A) - Fy^2{^) is used to estimate the 
uncertainties of the observation. However, d{A) is not the correct 
error estimator; although it does contain the errors due to statis- 
tical fluctuations, it is also sensitive to systematic effects such 



Table 1. Object sample and stellar properties (X-ray luminosity 
Lx , hardness H, spectral type, photosphere temperature and 
stellar age) published in Giidel et al. ( 2007) 



Name 




Hi 


Spect;_ 




Age"* 




[lO^o erg/s] 






[K] 


[My" 


04187-1-1927 


0.91 


0.57 


MO 


3850 




04303-1-2240 


4.99 


0.76 


M0.5 


3700 


0.5 


04385-f2550 


0.40 


2.26 


MO 


3850 


— 


AATau 


1.24 


1.98 


K7 


4060 


24 


BPTau 


1.36 


0.68 


K7 


4060 


1.9 


CITau 


0.19 


1.62 


K7 


4060 


2.2 


CoKu Tau/3 


5.83 


0.68 


Ml 


3705 


0.9 


CWTau 


2.84 


_ 


K3 


4730 


7.0 


CYTau 


0.13 


0.46 


M1.5 


3632 


1.5 


CZTau 


0.42 


0.43 


M3 


3415 


2.1 


DDTau 


0.09 


2.55 


M3 


3412 


4.5 


DKTau 


0.91 


0.75 


K7 


4060 


1.3 


DNTau 


1.15 


1.01 


MO 


3850 


1.1 


FMTau 


0.53 


1.34 


MO 


3850 


2.8 


FOTau 


0.06 


1.13 


M2 


3556 


1.5 


FQTau 


0.12 


0.36 


M3 


3416 


2.8 


FSTau 


3.21 


0.73 


MO 


3876 


3.1 


FVTau 


0.53 


0.50 


K5 


4395 


4.8 


FXTau 


0.50 


0.52 


Ml 


3720 


0.9 


FZTau 


0.64 


0.87 


MO 


3850 


1.1 


GHTau 


0.11 


0.98 


M1.5 


3631 


2.0 


GITau 


0.83 


0.90 


K7 


4060 


1.8 


GKTau 


1.47 


0.96 


K7 


4060 


1.2 


GNTau 


0.78 


1.18 


M2.5 


3488 


1.0 


GOTau 


0.25 


1.09 


MO 


3850 


3.8 


Haro6-13 


0.80 


1.99 


MO 


3800 


0.6 


Haro 6-28 


0.25 


0.75 


M2 


3556 


10 


HKTau 


0.08 


4.19 


M0.5 


3778 


1.8 


HOTau 


0.05 


0.08 


M0.5 


3778 


9.1 


HPTau 


2.54 


1.02 


K3 


4730 


6.9 


IQTau 


0.41 


1.17 


M0.5 


3778 


1.1 


ISTau 


0.66 


0.48 


K7 


3999 


4.2 


ITTau 


6.47 


2.05 


K2 


4900 


4.8 


MHO-3 


0.46 


1.26 


K7 


4060 


2 


RYTau 


5.50 


1.41 


Kl 


5080 


2.1 


UZ Tau/e 


0.89 


0.49 


Ml 


3705 


2.3 


V410 Anon 13 


0.01 


0.68 


M5.8 


3024 




V710Tau 


1.37 


0.76 


M0.5 


3778 


1.7 


V773 Tau 


9.46 


1.08 


K2 


4898 


6.4 


V807 Tau 


1.05 


0.52 


K7 


3999 


1.5 


V955 Tau 


1.62 


0.46 


K5 


4395 


6.7 


XZTau 


0.96 


1.05 


M2 


3561 


4.6 



" From Table 6 in lGudel et all d2007h 
* Hardness, see definition in Eq. {3) 
From Table 9 in Gudel et al.' ( ^OOj) 
From Table 10 in ,Gudel et ah ( .20"07h 



as the truncation of the flux by misalignment of one of the two 
observations. Therefore, we make a first order correction of the 
uncertainty by subtracting the averaged difference 



d(A) = d{A) - {d(A)} 



(1) 



where the average is taken over the wavelength range which is 
relevant for the fit procedure used later for the 10 /vm-silicate 
feature. This correction brings d{A) closer to the truly random 
noise contributions. However, d(A) can underestimate the purely 
statistical fluctuations (noise) by chance at various wavelengths. 
But assuming that the noise remains similar across a narrow win- 
dow in wavelength, we can obtain its statistical value from the 
fluctuations of d(A) itself. We use = 1 1 data points around 
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each wavelength bin Aj to derive the measurement uncertainties 
cr by computing a geometrical sum of their average and standard 
deviation of the systematic differences d{Ay. 



cr' 



(d(Aj) - {d(At))keAf 

N-l 



{{d{Ak)}keAf 



jeA 

with A = [i - 5, / + 5] 



(2) 



The X-ray data were obtained from the XMM-Newton 
Extended Survey of the Taurus Molecular Cloud (XEST) which 
consisted of 28 exposures in total, spread over the whole Taurus 
region. Parameters used for this work are the stellar X-ray lumi- 
nosity Lx and further, we defined the hardness H of the stellar 
X-ray emission as the ratio between the hard and soft luminos- 
ity components. This allows to measure the relative contribution 
of different energy bands to the X-ray radiation; the hard and 
soft bands comprise the 1-10 keV and the 0.3 - 1 keV ranges, 
respectively: 



H = 



Lx(l keV < £ < lOkeV) 
Lx(0.3 keV < £ < 1 keV) 



(3) 



As a few objects such as DH Tau show ed excessive emissio n 
during the observation (as described by ' Telleschi et al.l l2007h . 
their value for Lx does not correspond to the actual averaged lu- 
minosity which is the parameter of interest for the present study. 
Therefore, these objects have been discarded from our sample. 

The last colu mn of Table [T] lists the stellar age (from 
iGiidel et al.l2"007l) . These values were derived from T^g and L» as 
given in the literature and provide best-estimate values for which 
systematic uncertainties are difficult to provide. Although Giidel 
et al. indicated conservative age uncertainties (based on a variety 
of literature values for Teff and L*) of factors of 2-3, these are 
extreme values, and most a ges are - within t he framework of one 
set of evolutionary tracks (S iess et al.ll2000l in lGudel et al.ll2007l) 
- much more narrowly confined. As for t he original parameters 
used for age determination, Tf.fi and L», Whit e & GhezI (1200 Ih 
(used in the XEST stud y) estimate uncertain t ies in L* of 0.1- 
0.3 dex for their sample. iHartigan & Kenvoril (|2003|) (also used 
in the XEST study) estimate errors in L, from the disagreement 
of ages between components of binaries, amounting to typically 
0.1-0.2 dex. A similar uncertainty (0.11 dex) for L» has also 
been given bv lKenvon & HartmannI (IT995i) . A scatter of 0.2-0.3 
dex is furthermore found when comparing values from various 
authors. We thus conclude that most of our ages (primarily de- 
termined by L, and much less by T^ff) are accurate to within a 
factor of 1.5-2. 

The range of the stellar ages spreads from 0.5 to 10 Myr. 
Although our sample consists of Class II T Tauri stars only, this 
wide range indicates that the selected objects can be categorized 
into three evolutionary groups: Very young objects (<1 Myr) 
which are likely to be in a transitional phase from embedded 
to disk-only geometries; typical Class II objects with a pure 
disk geometry; and older objects (>5 Myr) which are more ma- 
ture Class II sources and are likely to have changed character- 
istics compared to their younger counterparts. These groups are 
not separated sharply given the age uncertainties. Also, we note 
that age is not the only parameter determining the evolutionary 
state of the disk. Therefore we continue to study the full sample 
and will investigate the object selection later in this paper (see 



3. The decomposition of the 10 jum silicate feature 

3.1. Modeling of the Emission Profile 

Our model describes the total observable flux with three com- 
ponents based on the two-layer temperature distribution (TLTD) 
method introduced by Juhasz et al. (2009): The emission from 
the stellar photosphere, a continuum emission from the opaque 
disk midplane and emission from the optically thin disk atmo- 
sphere consisting of thermal emission from dust grains of differ- 
ent mineralogy and temperature. The atmosphere is transparent 
with respect to the continuum emission of the disk midplane. 
Figure [U shows a sketch of the disk model. 




Fig. 1. Circumstellar disk model describing the observable total 
flux which is a superposition of the stellar light, F*, the contin- 
uum emission from the disk midplane, Fmid, and the emission 
from the optically thin disk atmosphere, Fatm- 



The total observable flux is given by 



^'v.tnf — f v.* + ^'v. 



(4) 



We describe the flux of the stellar photosphere approximately 
with the emission of a blackbody with a temperature J*: 



Fv,* - Ci,- 



2/1 



1 



ghv|kT^, _ I 



(5) 



C* is the normalization factor used for the later fit and By(T) 
defines the Planck function at temperature T. To model the flux 
from the midplane and the atmosphere of the disk, we set up two 
simple continuum emission profiles, one of which we use di- 
rectly for the midplane while the other is multiplied by dust grain 
absorption coefficients Kyj for modeling the flux from the atmo- 
sphere. These two continuum emission profiles are constructed 
as a superposition of blackbody spectra taking a radial temper- 
ature distribution T(r) and geometrical aspects of the disk into 
account, assuming axisymmetry. Hence, the fluxes can be writ- 
ten as 



Bv(T^id{r)) r dr 



By(TMr)) r dr 



(6) 



(7) 



where r is the radial distance to the central object, /?atm/mid,in/out 
the inner and outer radii of the disk midplane and atmosphere. 
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respectively (see Fig. [T]), C, are the normalization factors of the 
dust species / (in our case we use = 10, see Sect. |3.2| i and 
Co of the disk midplane emission, respectively. We assume the 
temperature distributions to follow a simple power law with 



(8) 
(9) 



7'mid('") 


( ' 


^mid,max 


\^mid,in 


ratm(^) 


( ' 


^atm,max 


\ ^atm,in 



which implies that all grains in the atmosphere or in the disk 
mid-plane follow the same temperature profile, regardless their 
size and chemical composition. This allows a substitution of r 
with T and Eqs. (|6]l and (|7]l can be rewritten as 



v,mid — 



Co 

^mid 



v,atm 



1=1 



^atm 



■r 



Bv(r) • r 



(2-?mid)/9mid 



dr 



(10) 



Bv(r) ■ r*^"''"""'^*'"" dr (ii) 



where C, are the new normalization factors used for the later fit. 
In this fitting approach, C*,Co, . . ., Ca-, ^^id, ^atm, T*, Tmid.max, 
rmid,min, T'atm.max and ratm.min are fitting parameters. As we are 
only interested in a wavelength range between ~7 yum and 14 yum 
(see Sect. |5.1| for a discussion), we can eliminate several of these 
parameters: 

- T-t, can be derived from the Uterature. We used the values 
summarized in Table [T] 

- We chose to set rmid.min and ratm.min to 10 K to account for 
the contribution of the cooler region of the disk which is ir- 
relevant for the mid-infrared regime. 

- Figure|2]shows the impact of Tmid.max on the shape of the flux 
function given by Eq. (fTOl i by keeping ^mid constant. It is ob- 



= -0.5 




Fig. 2. Normalized flux as expressed by Eq. (fTOl i for q - -0.5 
and Tmidmax = 250 K, 500 K, . . ., 2500 K. 



vious that Tmid.max has no significant influence on the shape 
of the flux function for values Tmid.max > 1000 K, which 
is true for various values of ^mid- Further, iD'Alessio et al.l 
([l998) showed that the inner disks of classical T Tauri stars 
(CTTS) reach temperatures around T^ub = 1800 - 2000 K 
due to the dust sublimation. Therefore, we set without loss 



of generality Tmid ,r 
tion we set T,,,^ 



1800 K. With the same argumenta- 



= 1800 K. 



- First trials of fitting simulated spectra showed that the fit 
is not sensitive to g'atm for any physically meaningful value 
(e.g., > ^atm ^ -2) as the influence on the shape of the flux 
function in Eq. ( fTTT l is fully dominated by the dust emission 
profile. We therefore set g'atm = ?mid = 

Consequently, the remaining fit parameters are C*, Co, . . ., CAf 
and q. For a given q and more than N + 2 data points, C*, Co, 
■ ■ Cn can be derived analytically according to a non-negative 
least-square fit (see e.g. lLawson & Hansonlll974 ) and further, a 
unique solution exists. Therefore, we can search efficiently for a 
global minimum ofx^(q). 

To derive the components C,- and q of the fit function and 
their errors for a given observation, we compute 100 spectra 
adding Monte-Carlo simulated, normally distributed noise to the 
original data; the noise distribution is based on the spectral er- 
rors of the original data. Each of these spectra is then fitted with a 
least-square method approach as we take the median values (the 
distributions tend to be very asymmetric) of the resulting param- 
eters to obtain the final parameters and their errors (Icr range of 
the frequency distribution). 

3.2. Emission profiles 

We use emission profiles analogous to 'Schegerer et al.' (l2006h 
and Bouwman et al. (2008): For the amorphous silicates we 
use profiles calculated for homogeneous, compact, and spheri- 
cal grains applying Mie theory. For the crystalline silicates we 
use emission profiles calculated for inhomogeneous s pheres ac- 
cording to the distribution of hollow spheres (DHS. iMin et al.l 
2005). We start from the complex refractive indices «, for sil- 
icate material /. The result is the dimensionless absorption ef- 
ficiency Qi, which is used for calculating the mass absorption 
coeflicient k^j - Qina^/(4/3na^p) where a is the particle ra- 
dius and p the material density. We fit the 10 //m silicate fea- 
ture with amorphous silicates with the stoichiometrics of olivine 
(MgFeSi04) and pyroxene (MgFe[Si03]2) and the crystalline 
silicates forsterite (Mg2Si04), enstatite (MgSiOs) and quartz 
(SiOa). Table |2] summarizes the dust species considered here. 



Table 2. Dust species used for the decomposition fit of the 1 jum 
feature fit. 



Name 


Stoicliiometry 


Structure" 


Reference 


Olivine 


MgFeSi04 


a 


Dorsciiner et al. 1995 


Pyroxene 


MgFe[SiO,]o 


a 


Dorschner et al. 1995 


Forsterite 


Mg2Si04 


c 


Servoin & Piriou 1973 


Enstatite 


MgSiOj 


c 


Jiiger et al. 1998"' 


Quartz 


SiO, 


c 


Seltzer & Kleinman 1960 



" a=amorphous, c=crystalline 



Fig 
(l2006l) 



Figure[3]shows a reproduction of Fig. 3 from i Schegerer et all 
where the mass absorption coefficients for different grain 
sizes and grain compositions are shown. 

In this work, we use only grains of sizes 0. 1 fim and 1 .26 pm 
as this is sufficient to fit the data reasonably well. Further, the 
fit with only two grain size s has been confi r med to be valid by 
iBouwman et~an (1200 lb and ISchegerer et al.1 (l2006h . The 10 /im 
silicate feature does not put constraints on silicate particles larger 
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E 
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MgFeSiO^ A 
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// 
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// 













10000 
8000 - 
6000 - 
4000 - 
2000 - 



Crystalline 




MgFe[SI0j]2 




8000 






6000 


// 




4000 


// ^ 




2000 



MgSiOj 




J 


V A 

" -■' " "1 ft 



SiO, 




8 • 





Fig. 3. Mass absorption coefficients of grains with radii 0. 1 yum 
(sofi dUne), 1 .26 //m (dashed line) and 2.5 jim (dotted lines). Plot 
after ISchegerer et"an (l2006l) . 



than about 5 /im in general. Therefore, discussing the crys- 
tallinity implies that we understand it as a fraction of sub- or 
micron sized grains. We do not fit features from PAH molecules. 
The latter produce emission lines at 7.7 yum, 8.6 fim, 11.2 fim 
and 12.8 //m with in the spectral range of our interest (see, e.g., 
iGeers et si. I I2006I) . Rather, we concentrate only on the silicates 
described above and perform the fit with 5 diff'erent silicates of 
2 diff'erent grain sizes each (consequently, = 10 different fit 
components). Therefore, we exclude wavelength regions in the 
spectra that clearly show PAH emission features . 

4. Results 

We list the fit results, i.e. the relative mass fractions of the miner- 
als as well as the derived values for ^disk, in Appendix|A] Table[3] 
lists values of the crystallinity F which is the sum of all relative 
mass fractions of the crystalline components regardless of their 
size. The reduced x^^^ values are listed in the third column, cal- 
culated from the averaged fit parameters. In Fig.|4]we present the 
spectrum and resulting fit functions for the example of AA Tau; 
the complete sample is shown in Appendix B. 

Most of the spectra are fitted reasonably well with a re- 
duced x^^^ in the range of 1-3. The fit of a few spectra such as 
IRAS 04303+2240, FV Tau, Haro 6-13, MHO-3, and RY Tau 
show systematic discrepancies between the data and the fit, re- 
sulting in a large x^^d- most of these cases, the data show a 
high signal-to-noise ratio and consequently, the error bars of the 
data are small. Our fit model is too incomplete to describe these 
object accurately. In the particular case of IRAS 04303+2240, 
the spectrum shows a large variety of emission features which 
we are not able to describe with the selected minerals. We did 
not intend to increase the complexity of the model to avoid in- 
cluding too many degree of freedom for the bulk of the sample 
with lower signal-to-noise ratios. We do not use these poorly fit- 
ted objects for our further studies. 



Fig. 4. IRS spectrum (black line) and fit (red line) including the 
corresponding uncertainties (gray lines) of AA Tau as an exam- 
ple of the full data sample. The continuum background used for 
the fit function is shown in blue. In the lower part of the fig- 
ure, the resulting Ax^, multiplied by the sign of the deviation, is 
shown. 



5. Discussion 

5. 1 . Disk model validity 

The fit method presented here for decomposing the 10 fj.m sil- 
icate feature takes advantage of a continuous temperature dis- 
tribution. Therefore, this m odel is more realistic t han previous 
methods using polynomials (;Bouwman et al. "20011), single tem- 
perature blackbodyjunctions (Meeus et al. 2003} ) or two tem per- 
ature fits (e.g. Ivan Boekel et al.ll2005l oriSargent et al."20 09|). For 
a syste matic comparison between the methods, see Juhasz et al.l 
(l2009l) . 

We conclude from the spectral fits that the continuous tem- 
perature distribution of the disk is a robust model for approxi- 
mating the continuum emission. We are able to adjust the back- 
ground by just one geometrical parameter, i.e. the slope q of the 
radial temperature distribution, which shows the strength and the 
simplicity of this method. Although previous studies pointed out 
that the exact shape of the background function has little effect 
on the relative composition of the dust mineralogy, the applica- 
tion of a more physical model is appropriate and avoids wrong 
conclusions about dust temperatures from a single blackbody ap- 
proach. Further, using the TLTD method to describe the silicate 
emission features of the disk atmosphere allows to model the 
co ntribution to t he 10 /im flux by grains of varying temperatures. 
As l juhasz et al.l S200% pointed out, this is crucial to reduce the 
systematic uncertainty of the decomposition. 

On the other hand, it has to be mentioned that the tempera- 
ture profile is assumed to be independent of the dust grain size 
and composition, which is an evident simplification, especially 
at the disk surface. Further, the TLTD method in its presented 
form does not account for radially dependent distributions of 
the individual dust species. As suggested bv lJuhasz et al.l (I2009h . 
to optimize the validity of the applied method, we confined the 
wavelength range to the 10 fj.m feature only and did not include 
longer wavelengths. It would be very interesting to extend the 
wavelength range to the full Spitzer IRS spectral coverage; a 
radially dependent distribution of the dust composition will be 
implemented in the TLTD method in a future study. 
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Table 3. Derived quantities for the decomposition analysis of 
the 10 fim silicate feature. The values in brackets correspond 
to the border values of the Icr range. The last two columns list 
the values of the crystallin ity for the cold (Fc ) and warm (Tw) 
components as derived bv iSargent et al.l (l2009l) . 



Name 




F r%i 


, 

A red 




rw''r%i 


04187+1927 


33.4 


(32.6, 33.6) 


2.7 






04303+2240 


33.6 


(32 8 34 4) 


14.3 






04385+2550 


1.7 


(1.5, 1.9) 


1.8 






AA Tau 


15.6 


(14.4, 16.9) 


1.3 


25±21 


5+5 


BP Tau 


8.8 


(8.1,9.1) 


0.7 


15+5 


18+6 


CI Tau 


15.3 


(14.4, 16.3) 


1.7 


8+13 


6+4 


CoKu Tau/3 


22.7 


(21.4, 24.2) 


1.8 


18±7 


11+4 


CW Tau 


20.7 


(20.1,26.4) 


2.3 


15+21 


6+9 


CY Tau 


37.0 


(32.7,43.1) 


1.0 


11+52 


30+27 


CZ Tau 


2.2 


(2.2, 2.3) 


0.8 






DD Tau 


17.6 


(16 6 18 7) 


2.2 


19+13 


20+10 


DK Tau 


17.3 


(16.7, 17.6) 


2.2 


56+14 


14+3 


DN Tau 


44.3 


(41.2,47.2) 


1.1 


7±17 


50+27 


FM Tau 


1.2 


(1.0, 1.3) 


1.2 


13+9 


4+5 


FO Tau 


24.7 


(19.1, 34.5) 


0.4 


14+7 


17+11 


FQ Tau 


50.7 


(50.1,52.3) 


1.1 


53±153 


14+8 


FS Tau 


1.1 


(1.0, 1.2) 


1.3 


9±10 


7+8 


FV Tau 


1.5 


(1.1,2.0) 


4.9 


6±7 


8+5 


FX Tau 


4.6 


(3.9, 5.4) 


1.9 


15±6 


14+4 


FZTau 


43.1 


(42.1,44.1) 


2.7 


49±28 


22+7 


GH Tau 


34.9 


G 1 8 39 0") 


1.4 


19±11 


12+8 


01 Tau 


10.3 


(10.0, 11.4) 


0.5 


13±7 


3+4 


OK Tau 


7.7 


(7.6,7.8) 


1.9 


26±9 


8+3 


GNTau 


17.2 


(16.4, 18.9) 


2.2 


64±19 


12+4 


GO Tau 


16.1 


(13.5, 19.6) 


1.1 


14±12 


9+6 


Haro 6-13 


1.4 


(1.4, 1.5) 


3.8 






Haro 6-28 


30.3 


(30.1, 31.3) 


1.0 


18±11 


45+13 


HK Tau 


40.8 


(37.3,44.5) 


1.6 


8±7 


24+8 


HO Tau 


5.0 


(4.4, 5.5) 


1.2 


17±23 


14+7 


HP Tau 


4.1 


(4.o! 4.4) 


1.6 


12±24 


4+4 


IQTau 


6.9 


(6.4,7.3) 


1.5 


7±31 


9+5 


IS Tau 


29.0 


(26.6,31.7) 


1.7 


100+45 


14+4 


IT Tau 


32.3 


(28.3,40.1) 


1.7 


8±11 


63+45 


MHO-3 


1.5 


(1.4,1.6) 


13.6 






RYTau 


7.0 


(6.7,7.5) 


6.6 






UZ Tau/e 


11.6 


(9.4,11.6) 


1.1 


5±10 


2+4 


V410 Anon 13 


59.1 


(54.6,61.2) 


1.6 


6±7 


43+13 


V710Tau 


16.2 


(15.3, 19.0) 


1.4 


37+21 


23+10 


V773 Tau 


22.4 


(22.4, 26.0) 


2.7 






V807 Tau 


32.1 


(28.6, 37.1) 


1.7 


15+7 


8+18 


V955 Tau 


36.5 


(33.8,38.7) 


1.9 


21+9 


25+8 


XZTau 


23.5 


(22.0, 25.4) 


1.8 


4+14 


8+12 



" From lSargent et alj l l2009l) 



Most obje cts of our sample have been studied likewise by 
ISarpent et al. l (2009) (SA09), using a two temperature decompo- 
sition fit (2T fit). Although this method appears to be less realis- 
tic, it probes different regions of the disk and is worth to compare 
with our results. The last two columns of Table[3]hst the accu- 
mulated values for all crystalline mass fractions for the warm 
and the cold fit components, respectively. The errors have been 
calculated using standard error propagation of the uncertainties 
given in Table 6 of SA09. The method used by SA09 to derive 
the uncertainties of the fitting parameters is dissimilar from our 
work and the values for the uncertainties as presented by SA09 
are very conservative. The consequence is that large errors ap- 
pear in the propagated values for the crystallinity as shown in 
Table[3] Therefore, we expect many objects to have overlaps be- 



tween the 1 sigma range of our F-values and the values derived 
by SA09. Indeed, out of 34 common objects, 23 show an over- 
lap with either the warm or the cold component. This number 
increases to 28 objects when considering our results with a 2 
sigma uncertainty. However, such an ad hoc comparison is not 
very valuable as a real comparison of the derived crystallinity is 
not possible: Basically, we should compare our results with the 
warm dust component of SA09 only as this corresponds more 
likely to the 10 //m feature. However, the largest fraction of the 
mass contribution for the two temperature fit is contained in the 
cold component (~ 100 - 200 K) which does not contribute to the 
10 /vm flux significantly. 

5.2. Correlations between crystallinity and X-ray emission 

We compare stellar high-energy properties with structural char- 
acteristics of the dust disk. For this purpose, we correlate the 
X-ray luminosity Lx and the product of Lx and Hardness H 
with the total crystalline mass fraction F. Using the full sam- 
ple of the remaining 37 objects shows that these quantities do 
not correlate. As discussed in Sect. |2] our sample consists of 
objects of a wide spread in age and therefore it is likely that 
different evolutionary stages are present even if we are study- 
ing class II objects only. As mentioned in Sect. |2l to increase 
the uniformity of the evolutionary epoch, we selected the ob- 
jects according their age and divided them into three groups; 
very young, intermediate and older objects. Unfortunately, our 
sample size decreases to 34 as we lose another three objects 
(IRAS 04187+1927, IRAS 04385+2550 and V410 Anon 13) for 
which no age determination could be found in the literature. 

We found a significant correlation between Lx and F when 
selecting the objects in the intermediate range (~1 Myr to 
~5 Myr). Figure |5] shows the correlation plots for the three 
groups while the age limits of 1.1 Myr and 4.5 Myr have been 
set to optimize the correlation performance with the number of 
objects used for the intermediate group (see below). We measure 
a correlation coefficient of -0.62 and a significance of 99.7% for 
a correlation. 

We calculate the regression lines Lx(F) according the 
ordinary least sq uare (OLS) bisector method described by 
llsobe etal.l(ll990h to treat the two variables symmetrically while 
we use the logarithmic value of Lx . The dashed lines in Fig. |5] 
correspond to lines with slopes adding and subtracting the re- 
gression slope unc ertainty as calculated according to Table 1 of 
llsobe etaI.I(fT990l) . 

We note here that measurement errors (or errors derived from 
measurements) will not be used in the derivation of regression 
lines. If measurement errors are a minor contribution to the scat- 
ter around a regression line, then the deviation of a point from 
the best-fit line is dominated by systematic processes not consid- 
ered here; using measurement errors as weights for such points 
will introduce arbitrary bias that is unrelated to the actual scatter 
of the points. This problem has been discussed in detail by Isobe 
et al. (1990) and Feigelson & Babu (1992). Is the scatter domi- 
nated by unknown, additional processes in our case? The uncer- 
tainties in Lx of our sample were discussed in Giidel et al. (2007) 
with the result that the errors resulting from the X-ray measure- 
ment process and the spectral fit procedures are usually smaller 
than variations due to intrinsic X-ray variability on time scales 
of hours to days. The uncertainty introduced by X-ray variabil- 
ity usually corresponds to a factor of ~2 between maximum and 
minimum values. As can be seen in our correlation plots dis- 
cussed here, such variations are smaller than the scatter around 
the regression lines. We conclude that other factors not consid- 
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Fig. 5. X-ray luminosity Lx vs. total crystalline mass fraction 
r for objects of an age within 1.1 - 4.5 Myr (top panel), for 
objects younger than 1 . 1 Myr (bottom left panel), and older than 
4.5 Myr (bottom right panel). The lines in the top panel represent 
the OLS bisector regression (solid) with the uncertainties for the 
slope (dashed). 



ered here dominate the scatter, and errors from the measurement 
or spectral fit procedures are inappropriate to use here. We there- 
fore use unweighted regression (Isobe et al. 1990). 

To define the limits of the stellar age for which the cor- 
relation works best, we varied the lower and upper limit and 
calculated the correlation coefficient of the intermediate group. 
Figure |6] shows a map of the correlation coefficient and the cor- 
relation significance as a function of the minimum and maxi- 
mum stellar age used for selecting the sample. We see that the 
minimum and maximum stellar ages for which the correlation is 
still present, is not sharply defined: It varies between 1-2 Myr 
for the minimum and 2-4.5 Myr for the maximum age, respec- 
tively. This is in agreement with the age uncertainties discussed 
in Sect. |2l Beside the goodness of the fit, we also considered 
the statistics in terms of the number of selected objects which is 
shown in the upper right panel of Fig. |6l Since we tried to opti- 
mize the selection in terms of correlation and statistical sample, 
we show in the lower right panel of the same figure the product 
of the correlation coefficient and the number of objects. We de- 
cided to use objects of age older than 1.1 Myr and younger than 
4.5 Myr and 20 objects satisfied this condition while 14 objects 
fall outside the borders. We emphasize that these age limits are 
somewhat arbitrary within the typical age uncertainties adopted 
for our stellar ages (see Sect. 2) and correspond to an optimum 
choice based on the limited number of objects in our sample. The 
selected age interval should be interpreted as essentially contain- 
ing CTTS of typical ages in Taurus (^1-5 Myr). We marked 
the final selection with a circle in Fig.|6l 



Since we selected the sample according the goodness of de- 
compositional fit of the 10 //m feature using < 3 as a cri- 
terium, we investigated how this cut effects the correlation. For 
this purpose we added all objects to the sample within the opti- 
mum age limits but higher ^^^^j, in particular we added MHO-3 
and RY Tau. The correlation coefficient of this enlarged sample 
was found to be -0.63 with a probability of 99.6%. This is very 
comparable to the original result and we conclude that the ap- 
plied threshold for x^^^ has no influence on the systematics of 
our study. 

Finally, we plot in Fig. [T] the product of Lx and H against 
r, using the same object selection as before. We measure a cor- 
relation coefficient of -0.66 with a significance of 99.9%. The 
product of Lx and H may represent the deposited energy of the 
X-rays in the disk. It is remarkable that the correlation for the 
product of Lx and the hardness correlate better with the crys- 
talline mass fraction than the X-ray luminosity alone. We inves- 
tigated the optimum age selection limits as before, and got for 
both the minimum and maximum age identical values as derived 
for the Lx vs. F correlation. 

We repeated this study with the values of the crystalline mass 
fraction derived by SA09 to discuss the dependencies of the 
correlation on the applied method. Therefore, we correlated the 
crystalline mass fractions Fc and Fw (see Sect. 15. H and Table [3]) 
with Lx ■ Since the reduced;^^^^^ in SA09 are systematically high, 
we had to increase the limit to x^^^ < 5 to calculate the correla- 
tion with a sufficient sample. With this approach we use a sam- 
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Fig. 6. Maps of the correlation coefficient (absolute value, upper 
left), the number of selected data points (upper right), the signif- 
icance (lower left) and the product of the correlation coefficient 
and the number of selected datapoint (lower right) for the Lx vs. 
F correlation taking only objects into consideration which have 
a stellar age in between the minimum and the maximum age as 
indicated by the two axes of the maps. The circles indicate our 
selection of the stellar age for which the correlation coefficient 
and the statistics has been optimized. The straight lines indicate 
where the minimum and maximum age limitations are identical. 
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Fig. 7. Product of the X-ray luminosity Lx and the X-ray hard- 
ness H vs. total crystalline mass fraction F for objects of an 
age within 1.1 - 4.5 Myr (top panel), for objects younger than 
1 . 1 Myr (bottom left panel), and older than 4.5 Myr (bottom right 
panel). 



pie of 20 common objects of which 14 fall into the age range 
of 1.1-4.5 Myr. We found that r,,, correlates weakly with Lx 
(corr=-0.54, P=95.6 %) and strongly with H ■ Lx (corr=-0.65, 
P=9S.9 %). Fc does not correlate with any of these parameters 
at all. Further, repeating the investigation of the age range selec- 
tion showed a comparable picture to the valid range (1.1- 1.5 Myr 
for the lower limit and 1.5-6 Myr for the upper limit). 

This result confirms three aspects of our study: First, the cor- 
relation seems to be independent of the applied fitting method 
and is therefore not biased by the systematics of the decomposi- 
tion approach. Second, our TLTD fit refers to the 10 fim feature, 
for which we found a correlation between crystallinity and Lx 
. In the 2T fit (as used by SA09), the 10 fj.m feature is emitted 
by the warm component, and a corresponding correlation is in- 
deed found. Our focus on the 10 fim silicate feature is therefore 
justified a posteriori as the cold component from the 2T fit does 
not show a correlation. And third, we see in both data sets that 
the correlation for the combined product of hardness and X-ray 
luminosity is better than for the X-ray luminosity alone. 

5.3. Is the dust amorphized by the stellar wind? 

It is a remarkable result that an anticorrelation between the X- 
ray emission and the crystalline mass fraction is found; so far, 
no correlation between properties of the central object and the 
crystalline structure of the dust disk has been reported from ob- 
servations (e.g., Sicilia-Aguilar et al. 2007). 

However, the observed anticorrelation demands an indirect 
explanation because X-rays of the observed energies carry too 



little momentum to damage the crystalline structure of the dust. 
Their energy will rather be absorbed by the electrons, which tem- 
porarily leads to ionization of lattice atoms and finally to the pro- 
duction of phonons and a resulting increase of the temperature 
within the grain. 

Although the processes of X-ray emission are still contro- 
versial for T Tauri objects, it seems clear that they mostly orig- 
inate in magnetic coronae analogous to the solar corona, al- 
though stellar coronae are much more X-ray luminous and hot- 
ter. Therefore, we assume that the X-ray emission of our sample 
of T Tauri stars is related to the high energy processes in the stel- 
lar corona that also lead to a solar-like wind composed of ions, 
electrons and neutrons. We show in the following that the flux 
of the solar wind and the particle energies would be sufficient 
to amorphize dust grains efficiently in a protoplanetary disk at a 
radius of 1 AU and that this process induces significant changes 
in the emission profile of the dust grains. Although we may ob- 
serve the silicates at larger distances from their central object, 
the choice of 1 AU is based on the availability of data from the 
solar wind and allows for an extrapolation to larger distances. We 
use a simple disk model where homogeneous dust grains at the 
disk surface are directly exposed to the stellar wind and where 
recrystallization effects are neglected. 

5.3.1. Stellar wind properties 

We first derive the particle fluxes and energies from data from 
the solar wind as this is the only object for which direct and reli- 
able measurements of these parameters are available. We know 
from the enrichment of spallation products in meteorites that the 
young Sun had a proton flux several orders of magnitude higher 
than at present (CafFe et al. 1987). Unfortunately, the energies 
which lead to nuclear spallation are in the order of £ > 10 MeV 
while the energies o f interest for amo rphization are £ < 1 MeV 
(see below and, e.g.. |Jager et al.ll2003l) . The processes producing 
these energies might be dififerent and consequently it may not be 
correct to scale particle fluxes with the well measured soft X-ray 
flux for younger stars. Therefore, we proceed with present solar 
particle fluxes and energies and use its amorphization potential 
as a lower limit. 

We used the database from OMNIWeb jKing & Papitashvilil 
2005i) to determine the present proton density of the solar wind 
plasma, its velocities and proton to He-ion ratio near the Earth. 
We averaged the full dataset over the years 1963 until 2008 
and obtained a density of ~7.2 protons/cm^*, a wind speed of 
~450 km/s and a He-ion to proton ratio of ~4.4 %. Uncertainties 
of these quantities are up to 30 %. This leads to a mean proton 
flux of 3.2-10** protons/(cm^-s) at energies around 1 keV and to 
a He-ion flux of 1.4-10^ ions/(cm^-s) at energies around 4 keV. 



5.3.2. Calculation of dust amorphization 

We use these values to compare them with the ion dose required 
to amorphize a dust grain. For this purpose we calculate the num- 
ber of displacements of lattice atoms per incident ion using the 
SRIM-2008 software (Ziegler et al. 2008). SRIM is a collection 
of software packages which calculate many features of the trans- 
port of ions in matter such as ion stopping, range and straggling 
distributions in multilayer targets of any material. Further, it al- 
lows the calculation of ion implantation including damages to 
solid targets by atom displacement, sputtering and transmission 
in mixed gas/solid targets. Table |4] summarizes the calculation 
performed for two minerals, pyroxene and enstatite, respectively. 
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irradiated by protons and He-ions at various energies. With the 
SRIM software we determined the penetration depth (range) of 
the ion and the number of displaced lattice atoms (due to elastic 
scattering on the atom's nuclei) per incident ion. 

The mineral densities have been taken from the webmin- 
eral databas^D. We determine the required ion dose for complete 
amorphization of the dust grain material down to the penetration 
depth of the incident ion by deriving the column density of the 
lattice atoms. Further, we require a quantity for how many dis- 
placements per lattice atom (dpa) are required to observe amor- 
phization. We set this number to be 2.5 based on the mineral ir- 
radiation experime nts for enstat ite dJager e t al. 200 3|), forsterite 
( Bringa et al.l2007l) and olivine ( Demvk et al..2001;iCarrez et al.l 
I2002h . However, this quantity remains inaccurate as most of the 
above cited studies determined doses for He-irradiation to fully 
amorphize t he mi nerals varying within an order of magnitude. 
iJager et al] ( l2003l) presented a dpa-value of 2.5 but the lowest 
He-ion energy used was 50 keV. For protons, we assume the 
same dpa ratio although this is hypothetical. 

Table |4] further shows that A^oisp/ion increases with increas- 
ing energy but the required dose for amorphization increases as 
well. This can be explained from the definition of the required 
dose. With higher energies, the ions penetrate deeper into the 
dust grain and consequently a thicker layer will be amorphized. 
This requires a larger total number of displacements as more lat- 
tice atoms are impacted. 

We can therefore conclude from this calculation that the 
He-ions of the present solar wind would amorphize the top 
~30 nm of the irradiated face of a pyroxene dust grain and the re- 
quired dose of ~2.3- 10"' ions/cm^ would be accumulated within 
~50 years. To allow the full surface to be amorphized, the dust 
grain has to rotate with respect to the ion beam to allow for an 
isotropic irradiation, and the timescale has to be increased by a 
factor of a few. Obviously this indicates a very efficient mecha- 
nism for amorphization and might be even more efficient as we 
have ignored the irradiation by other ions. 

On the other hand, this calculation assumes that all dust 
grains are fully exposed to the radiation of the central object. 
This is obviously not true for the majority of the observable 
grains in the disk atmosphere. The ion flux gets extinct analogue 
to optical light due to absorption by dust grains (self-shielding) 
and therefore, the amorphization timescale depends strongly on 
the location of the grain within the disk and the presence of ver- 
tical mixing of the disk material. Therefore, the 50 years of irra- 
diation time represent the timescale for amorphization, based on 
the full flux. This number will increase exponentially (or faster) 
the deeper the grains are located in the disk. Consequently, our 
observations probe various timescales for dust amorphization; 
conclusions depending on whether the presented mechanisms 
are too efficient or too inefficient cannot be made at the present 
time. Knowing the mean free path length along the particle tra- 
jectories would enable more quantitative conclusions about the 
amorphization timescale and potential with respect to other dust 
processing mechanisms. Particles of higher energies could also 
be considered by studying multiple scattering of ions with dust 
grains. This would allow an amorphization of dust deeper in the 
disk. For this purpose, dedicated modeling of the disk geometry 
and grain size distribution is required. This will be addressed in 
future studies. 

For now, we can conclude that the dust amorphization by 
stellar wind ions at the disk surface is sufficiently efficient so 



that this mechanism might dominate the dust processing at the 
disk surface layer. 

5.3.3. Optical properties and grain size dependency 

Since the dust grains are amorphized only at the surface layer, 
the impact on the optical emissivity has to be investigated to 
demonstrate the observability of the amorphization of circum- 
stellar dust by the low energetic stellar wind. 

As an example, we calculate the optical mass absorption co- 
efficients of crystalline pyroxene, coated with a 30 nm thick 
layer of amorphou s pyrox ene using Mie theory according to 
iBohren & Huffmanl (Il983h for coated spheres. The c alculation 
code in the a ppendix B of iBohren & Huffmanl d 19831) has been 
translated by iMatzled (12002') into a Matlab script which we 
have u sed for this study. The optical constants from lJager et al.l 
were used for crystalline pyrox ene while parameters fo r 
the amorphous pyroxene are taken from lDorschner et"aLl ( Il995h . 

Figure [8] shows the resulting mass absorption coefficient 
within the wavelength range of interest for a dust grain of 0. 1 jum 
and 1.26/im in radius, respectively. Further, this plot indicates 



MgFerS'IOjlj ■ ' ■ ■ ■ ' ■ ■ 


r = 0, 1 /im 




r = 28.1% 










www.webmineral.com 



6 8 10 '2 1-1 

A [Mm] 

Fig. 8. The thick, solid line corresponds to the mass absorption 
coefficient of a 2-layer spherical dust grain with pyroxene sto- 
ichiometry and a radius of 0.1 ;um (left) and 1.26 /im (right). 
The outer layer is amorphous with a thickness of 30 nm in both 
cases and the core is crystalline. The dotted line corresponds to a 
pure crystalline and the dashed line to a pure amorphous grain of 
the same size and chemical composition. The red, thin solid line 
shows a decomposition fit of the 2-layer absorption coefficient 
with a purely amorphous and purely crystalline component and 
a crystalline mass fraction of r=28.1% and r=93.5% is mea- 
sured, respectively. Due to the high crystalline mass fraction for 
the larger grain, the dotted, red and thick lines overlap in the 
right panel. 



the two extremes of a purely crystalline grain and a purely amor- 
phous pyroxene grain. Finally, we calculate a linear combination 
of the two extremes to fit the 2-layer case and get a contribution 
ratio for the crystalline dust of r=28.1% for the 0.1 jum grain 
and r=93.5% for the 1.26 fim grain, where F corresponds to 
the same definition as used for the decomposition fits. This re- 
sult reflects the mass ratio used for the coated sphere calcula- 
tion (29.1% crystalline for the 0.1 jum grain and 93.0% for the 
1 .26 jum grain) with relatively high accuracy. Consequently, even 
if the dust grains are amorphized in the surface layer only, the 
spectrally derived value for amorphization corresponds to the 
mass ratio of crystalline and amorphous silicates and conclu- 
sions about dust amorphization are valid even for low-energy 
ions. 
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Table 4. Ion-penetration depth (range), number of displaced lattice atoms per ion A^oisp./ion, and the minimum required dose for 
full amorphization of the upper layer (of thickness equal to the ion penetration depth) of pyroxene and enstatite by proton and 
He-ion irradi ation at variou s energies. For this calculation, mineral data from the webmineral data base and the calculation software 
SRIM-2008 (iZiegler et al.l l2008) have been used. For the calculation of the dose, a minimal displacement per lattice atom (dpa) of 
2.5 have been assumed to observe amorphization. 
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However, the provided calculation assumes a purely crys- 
talline dust grain. We expect to see a mixture of amorphous and 
crystalline dust regardless of the stellar wind and its impact on 
the dust structure. Consequently, the model described here deliv- 
ers only an additional component for the amorphous dust content 
and values for F are expected to be lower than these limits. 

Further, the results of the decompositional fits of our sam- 
ple (see Table [All suggest that a significant amount of the dust 
is present in larger grains and hence, the discussed mechanism 
might not be applicable. However, we repeated this study con- 
sidering only small grains and no correlation between Lx and F 
was found. It is very likely that higher energetic ions are im- 
portant and more abundant than in the present solar system. 
These ions penetrate deeper into the dust grain and consequently, 
larger grains are amorphized. Therefore, the approach of using a 
monochromatic ion beam with solar wind properties is too ele- 
mentary and has to be extended for more quantitative studies. 

6. Conclusions 

The previous calculations show that the amorphization of the 
surface layer of the protoplanetary dust by the stellar wind is a 
viable scenario. We have seen that a stellar wind with solar char- 
acteristics allows a very fast amorphization of the upper layer of 
dust grains compared to the longer timescales of dust processing 
during the star and disk formation. For sub-micron grains, this 
process leads to an observable structural change of the grain. 

However, these processes are inefficient for micron-sized 
grains or larger bodies as the penetration depth of the low- 
energetic ions are too small to degrade the dust sufficiently. 

On the other hand, as young stars are generally more active, 
we can assume that particle fluxes and energies are generally 
higher than in the present solar wind and consequently, the de- 
scribed processes may be even more efficient. Further, the solar 
irradiation is not only composed of the solar wind but a great va- 



riety of additional particles at various energies and fluxes. Since 
our knowledge of these particle fluxes (E %\Q MeV) in the early 
time of the solar history is very poor, and since particle irradi- 
ation in other systems are unmeasurable, the assessment of the 
amorphization efficiency remains very inaccurate and difficult. 
Other irradiation sources than the star, such as shocks from stel- 
lar winds in the protoplanetary gas disk or from accretion or jets 
confuse the overall picture further Recrystallizati on due to ther- 
mal annealing of warmer dust grains (see, e.g., iDiouadi et al.l 
2005) may inhibit the amorphization process and needs to be 
taken into account for quantitative models. 

Nonetheless, we observe a correlation for 20 objects between 
the X-ray luminosity as well as the X-ray luminosity multiplied 
with the X-ray hardness and the crystalline mass fraction of the 
atmospheric dust of the inner protoplanetary disk for objects 
within an age range of approx. 1 to 4.5 Myr We interpret this 
as an indicator of high-energy processes in the central object. 
We therefore postulate degenerative processes for the crystalline 
structure of the dust by ionic irradiation. Although we cannot 
observe the stellar wind directly, its flux and/or speed might be 
related to the X-ray luminosity, as shown for low X-ray fluxes 
by W ood et a l. (2005) by comparing stellar mass loss rates with 
X-ray surface fluxes. 

Self-shielding effects of the dust disk play an important role 
for the time scale and the overall potential of amorphization. 
Consequently, the disk geometry and the dust density distribu- 
tion have to be taken into account for studying evolutionary ef- 
fects. Further, self-shielding could explain why the correlation 
diminishes if we include objects younger than ~1 Myr. These 
objects might have well mixed disks even in the upper disk atmo- 
sphere where dust settling has poorly progressed; consequently, 
the self-shielding is more relevant. 

It is less clear why the correlation worsens if we include 
objects older than ~4.5 Myr It could be related to a statistical 
problem as our sample does not include many objects at such 
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ages. Further physical effects may also play a role on longer time 
scales. Possibilities include radial dust mixing or the reproduc- 
tion of crystalline dust material in the inner disk region. 

It would be interesting to verify if this correlation is present 
in other star forming regions and for a larger sample of objects. 
The extension of this study to Herbig Ae/Be systems could pro- 
vide different aspects of this process due to the shorter evolu- 
tionary time scales of the central object. 
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Tkble A.l. Resulting fit parameters for the decomposition analysis of the 10 fxm silicate feature. Dashes indicate that the lower boundary has been reached and the component 
proportion has been set to 0. 
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Appendix B: Spectra of the complete Data Sample 



Table B.l. IRS spectra (black line) and fit (red line) including the corresponding 
uncertainties (gray lines) of tlie full dataset. The continuum background used for 
the fit function is shown in blue. In the lower part of the figures, the resulting 
multiplied by the sign of the deviation, are shown. 
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